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Abstract 

A natural procedure is presented to measure the angle 7 from the decay — > 
. It is based in the Dalitz plot fitting analysis. Neither amplitudes nor 
strong phases have to be known a priori. We present simulations of this decay com- 
puting both statistical and theoretical uncertainties and analyze the experimental 
feasibility. We found that 7 could be measured with a combined error of the order 
of 20° with 90% of CL after about a couple of years of running of the first generation 
of B factories. 
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The study of hadronic decays in the B system seems to be a powerful tool for the 
understanding of CP violation. To check Standard Model (SM) predictions it is par- 
ticularly important to measure the three Cabibbo-Kobayashi-Maskawa (CKM) angles 
a = arg(-V td V t * b /V ud V: b ),P ee arg(-V cd V; b /V td V t * b ) and 7 = arg{-V ud V* b /V cd V* b ) . Only 
(3 is expected to be clearly measured from the gold plated B Q — ► J/ipK s decay which is 
almost free from theoretical uncertainties |1J and benefits from the large (10 -3 ) branching 
ratio. The extraction of the two other angles requires the measurement of decays with 
branching ratios of the order of 10~ 5 or less. Many interesting methods using various 
decays have been proposed so far in the literature P, ||, U but the matter is still open. 

More precisely, the angle 7 seems to be hard to measure. The method presented in 
Ref. @ provides a theoretically clean procedure to extract 7 combining decays with D° 
in the final state. Unfortunately, both the original method and clever extensions [Q] of it 
demand a large statistics. As a result, one expects to need about ten years of data 
taking in the first generation of B factories to attain a reasonable error — at least 15° 
- in the measurement of 7. It is then interesting to look for other methods that could 
provide a constraint for the value of 7 in much less time — their eventual theoretical 
errors should be as well estimated as possible 0. 

In this letter we present a direct and simple method that could provide a nice first 
measurement of the angle 7 after about a couple of years of data taking in the first gen- 
eration of B factories. We use the decay — > tt ± 7i + ii^ where the necessary interference 
is given by the intermediate resonant channel Xco^. This channel has been first pointed 
out in Ref. J7|; nevertheless, in that reference the method used to extract the angle was 
very model dependent and demanded large statistics. Here, we present a totally different 
approach: we show the viability of performing a full Dalitz plot analysis of this decay. It 
can provide a direct measurement of the angle 7 free from model dependencies. 

Other methods existing in the literature to measure 7, independently of the considered 
channel, are based in the measurement of branching ratios and asymmetries 0. The 
relationship between these measured numbers and the angle 7 is not direct. Moreover, 
these methods generally present discrete ambiguities. 

The main feature of our method is that it exploits the fact that in three body decays 
one can have a direct measurement of the amplitude of a decay — instead of branching 
ratios, that is amplitude squared. This means that one can have a direct experimental 
access to the phase of a given decay. This fact has already been used|| in connection 
with CP violation, in a quite different context. The method presented in this letter can 
eventually be also used to extract CP violating angles from other three body decays of 
charged B's. 

Let us present our ideas using the channel B + — > 7t + 7t + tt~. Many intermediate chan- 
nels contribute. Indeed, resonant channels — p ^ + , /o 7r+ ) Xco 7r+ ; e ^ c — together with 
the direct non resonant decay produce the same experimentally detected final state. This 
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final state is thus the product of the interference of all these intermediate states. 

The fact that in three body decays one can measure differential widths — usually 
displayed in a Dalitz plot (DP) — allows a clean separation of these partial channels. 
The distribution of measured events in the plot can be fitted using appropriate fitting 
functions. 

The fitting technique has proven to be very successful in describing, for example, three 
body decays of D mesons||, even with only about one hundred of reconstructed events flO| . 

In order to do so, one considers a fitting function including one term for each interme- 
diate channel contributing to the final state. For example, for the decay B + — > 7r + 7r + 7r~ 
the fitting function should be 

F B +->ir+K+n- (ml, ml) = \Y,iaie tei Fi(m\,ml)\ 2 , (1) 

where m\ = (p 7r + + P-k-) 2 and m\ = (p 7r + + Pn-) 2 are the usual Dalitz plot invariant 
variables, Fi are the amplitudes corresponding to each partial channel and dj and 9i are 
unknown real parameters that will emerge from the fit. The sum is performed over all the 
intermediate resonances as well as the non-resonant decay. For the resonant channels, the 
function F± is very well known: it is simply the usual Breit-Wigner |TTJ times an angular 



function according to the spin of the resonance. The non-resonant decay amplitude is 



discussed in ||12|| ; we will be back to it later. 

The Dalitz plot maximum-likelihood technique uses the function of Eq. (j^) to fit the 
measured differential width distribution dYjdm\dm\ of the total decay. The output is 
the amplitude fractions a, and phases $i of each partial decay. In other words, it brings a 
direct measurement of all the phases. 

In general, these phases can be written as 9i — 8i + <fii, where 5i is a CP conserving 
and fa is a CP violating phase, respectively. Obviously, in this way it is not possible to 
separate 8i from fa because only their sum is measured. 

Nevertheless, now consider the CP conjugated decay B~ — > tt~7t~tt + . The phase of 
each partial amplitude changes to 9i = 5i — fa. 

If one then applies the fitting procedure to the DP's corresponding to both B + and 
B~ decays, one gets a direct measurement of the CP violating phase, 

fa = (6i - 9i)/2 . (2) 

This procedure does not require to make any assumption about FSI as other methods 
to measure CP violation demand[|TJ. In B meson decays, FSI are usually assumed to 
be small [|13|], but this is not necessarily correct |T^, [15]]. Here, strong phases Si do not 



have to be known a priori. Moreover, using this procedure one can also obtain a direct 
measurement of the strong phases 5i = (9i + 9i)/2. This by-product of our method could 
be an interesting input to other methods. 
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For this procedure to apply one needs at least two intermediate channels with different 
CP violating phases. Indeed, if all the intermediate channels have the same CP violating 
phase it would factors out; in other words, there would be no interference to peak CP 
violation. 

In the decay — > 7r ± 7r + 7r _ , the Xco^ partial channel produces the necessary inter- 
ference to extract the angle 7. This channel is driven by the CKM coefficients Vb c V* d and 
it thus has no CP violating phase. On the other side, the direct non-resonant contribu- 
tion as well as the other resonant channels — p ^, fo^, etc — proceed via the CKM 
coefficients V u bV* d and their amplitudes thus contain the weak phase 7. 

Unfortunately, all the partial channels but XcqK contain also penguin diagrams which 
are driven by another CP violating angle, (3. These penguin contributions are expected 
to be small but not negligible [I], [L6[. 

In the following, we will present our simulations of the B ± decays. As a first step, 
we will not include penguin contributions. The corrections due to their inclusion will be 
studied later in this letter. They are the only source of theoretical uncertainties within 
the method presented here. 

The experimental simulation consisted in the following. First, we have generated a 
sample of B + — > 7r + 7r + 7r~ events using Monte Carlo technique. The dynamics was given 
by the function T of Eq. ([J), with a given set of input parameters and Then, we 
fitted the generated distribution of events in the DP using the maximum-likelihood fitting 
technique and MINUIT package JUjj . The fitting function was T but now a» and 9i are the 
floating parameters which are obtained from the fit. These two steps were then repeated 
for the CP conjugated decay B~ — > 7r~7r~7r + . 

In fact, we have considered a number of sets of input parameters Oj and $i correspond- 
ing to various possible scenarios for the unknown quantities involving this decay: the 
relative weight of each partial channel, their relative strong phases, and the angle 7. 

In Tables [l] and |2] we show the result of one of our simulations of the decay. It describes 
a probable scenario according to our present knowledge: 

1) BR(B+ -> Xc ^ + ) ~ 5 x 1O~ 5 0; BR(x c0 -> tt+tt-) ~ 0.8% @; BR(5+ -> 



TC + ic + ir) NR (non resonant) ~ HT 5 [|IJ; BR(5+ -»■ p°ir + ) ~ 8 x 1(T 6 |2§; BR(p° -> 
7T+7T-) ~100%. We have assumed BR(B + -> f ir + ) ~BR( J B+ p°n + ). From the square 
roots of the numbers above, one simply gets the coefficients ai of the column "input" in 
Tables | and @. 

2) The CP violating phase 02 = 03 = 04 = 7 for B + has been chosen as[^l|] 65°. The 



unknown strong FSI phases 8{ have been arbitrarily taken as 5°, 15° and —10° for NR, 
foir^ 1 and p ^, respectively. With these numbers one gets the values Q{ of the column 
"input" in Tables [I] and 0. 

The NR distribution has been considered flat. We fixed the Xco 71 " parameters a± = 1.0 
to have an overall normalization and Q\ = to fix our phase definition. 
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decay 


input 


200 events 


500 events 


1000 events 






1.0 


fixed 


fixed 


fixed 




Oi 


0° 


fixed 


fixed 


fixed 


NR 


a 2 


4.0 


3.1 ±0.7 


4.1 ±0.6 


3.8 ±0.5 




62 


70° 


(64 ± 25)° 


(71 ± 17)° 


(66 ± 11)° 


/OVT+ 


a 3 


2.5 


2.5 ±0.6 


2.8 ±0.4 


2.5 ±0.3 




O3 


80° 


(99 ± 28)° 


(88 ± 18)° 


(75 ± 12)° 




a 4 


3.0 


2.1 ±0.6 


3.4 ±0.5 


3.0 ±0.4 






55° 


(74 ± 26)° 


(39 ± 18)° 


(50 ± 12)° 



Table 1: Fitting results for B + Monte Carlo sample. 



decay 


input 


200 events 


500 events 


1000 events 






1.0 


fixed 


fixed 


fixed 






0° 


fixed 


fixed 


fixed 


NR 


a 2 


4.0 


4.2 ± 1.0 


3.6 ±0.6 


4.2 ±0.4 






-60° 


(-72 ±28)° 


(-64 ± 16)° 


(-66 ± 10)° 






2.5 


2.6 ±0.7 


2.4 ±0.6 


2.4 ±0.4 




h 


-50° 


(-48 ±29)° 


(-55 ±18)° 


(-48 ± 12)° 


p°7T- 


0,4 


3.0 


3.0 ±0.7 


2.8 ±0.7 


3.3 ±0.5 




h 


-75° 


(-81 ±27)° 


(-81 ± 16)° 


(-83 ±11)° 



Table 2: Fitting results for B Monte Carlo sample. 



We show in Table [l] the result of the simulation for the B + decay for three different 
numbers of generated events (200, 500 and 1000). In Table |2] we present the same results 
for the CP conjugated decay. 

One then uses Eq. (0) for the three CP phase changing channels NR, f 7i and pir. 
One gets 



7 = 


68 ±19 , 


7 = 


68 ± 12 , 


7 = 


66 ±7 


(NR) 


7 = 


73 ±20 , 


7 = 


71 ± 12 , 


7 = 


62 ±8 




7 = 


77 ±19 , 


7 = 


= 60 ± 12 , 


7 = 


= 66 ±8 





for 200, 500 and 1000 generated events, respectively. The errors in Eq. @ have been 
obtained summing in quadrature the independent errors of B + and B~ fits. 
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These results correspond to the scenario shown in Tables I and II. Nevertheless, as 
this scenario is based in particular assumptions we have performed a systematic study of 
these results, allowing a large variety of other scenarios. 

First, we have varied the BR of the partial channels - i.e., the square of the input 
coefficients a* - by as much as a factor of 5 and we got acceptable fits with similar errors. 
Second, we have tried other values of the CP conserving phases Si and we got the same 
accuracy for any value of the phases, even when they were all set to zero. Third, we have 
tried many different values of 7 between and In and we have always found the same 
accuracy in the results. Finally, we have made simulations releasing the shape of the 
function F2 describing the non-resonant channel [E| and found no important variations 
in the errors of Eq. (^) . 

We are then confident that in any acceptable scenario for this decay, the error to extract 
the angle 7 would be similar. This procedure thus brings a simple way of predicting the 
error in the measurement of 7; it would only depend on the number of reconstructed 
events — this is not the case for many other methods|5|]. 

It is worth mentioning another important point of our simulations. The method has 
no discrete ambiguities; accordingly, we always get only one value of 7 from the fit. 

Let us now discuss the real scenario including penguins, studying by how much their 
inclusion modifies our previous results. For example, in the fo^ + channel the measured 
quantity a 3 e ld3 is in fact 

a 3 e l63 = Te i(<5r+7) + Pe^ Sp ^ , (4) 

where Te l ^ T+ ^ is the tree contribution and Pe l( - Sp ~^ is the penguin one]l[; St and Sp are 
the strong phases. 

A pictorial representation of Eq. (|j) is shown in Figure 1. It shows that when measur- 
ing the angle #3 we are missing the actual tree phase by an angle e + . The same argument 
holds for B~, leading to an angle e_. As a result, Eq. (H) becomes 

03 = (^3-^3)/2 = 7 + (e + + e-)/2. (5) 

Thus, e = (e + + tJ)/2 is the theoretical error of our method. Figure 1 shows that the 
worst case corresponds to the configuration when the tree and the penguin contributions 
are orthogonal in the complex plane. As a consequence, we have 

|e±| < arctan(P/T) . (6) 

The actual value of the ratio P/T is not known at present. An estimate was obtained 
for the decay B — > nn; P/T ~ 0.2 [[16]. In our case, one expects the ratio P/T to be of 
the same order. Assuming P/T = 0.2, the uncertainty on 7 extraction due to penguin 
contribution would be at most ~ 11°. Anyway, Eq. (H) shows that as long as P/T remains 
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Figure 1: Tree, penguin, and the measured contributions. 



not very large, the penguin pollution does not invalidate the method; for example even 
for the improbable value P/T ~ 0.5, e < 26°. 

The inclusion of final state rescattering does not spoil our analysis of the error. The 
B + — > /o 7r+ decay proceeds through a unique isospin amplitude and thusfl5[] Eq. (£|) 
remains unchanged. For partial decays with more than one isospin amplitude — as B + — > 
p°7r + for example — the form of Eq. @ does not change either, but the interpretation of 
the two terms in this equation does. Indeed, even with rescattering one still has two types 
of diagrams: the first type having weak phase 7, the second having weak phase —j3. But 
now, for example, the term we call T in Eq. ([|) would be a more subtle combination of 
tree, color-suppressed and annihilation quark diagrams jn|. Thus, the method also applies 
to intermediate channels with more than one isospin amplitude; nevertheless, a complete 
isospin analysis would be required to stablish the theoretical uncertainty — note however 
that this is only necessary if rescattering effects are found to be large. 

Let us now study the experimental feasibility of this method. Using Eqs. (Q) and 
one can immediately get the error in the extraction of 7 according to the number 
of reconstructed events. For example, with 1000 reconstructed events and P/T = 0.2, 
this method would give 7 with a statistical plus theoretical error of 22° with 90% of CL. 
Moreover, combining the measurement of 7 using all the intermediate channels NR, pir 
and /o 71 " would certainly decrease the error of this procedure. 

Much less statistics is needed to simply detect CP violation in this decay. One observes 
CP violation when, e.g., 62 — 2 , is different from zero. Assuming 7 ~ 65° then from Eq. 
(H) one concludes that only 200 events are needed in order to detect a 3cr CP violation 
effect. 

Three main features allows us to be optimistic about the possibility of doing this 
analysis in a short period of time, either in BaBar, KEK or CLEOIII. First, tagging is 
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not required because one needs only charged B's. Second, the three detected particles 
are charged, thus we expect the efficiency to be high and the background to be not very 
large. Third, the method itself demands small statistics. For the B + — ► 7r + 7r + 7r~ decay, 
assuming for example a total BR of 3 x 1CT 5 and a reconstruction efficiency of 60% |J one 
would expect to need about 2 years of running of BaBar to reconstruct 1000 events. Of 
course, a full experimental simulation of this decay is required to have definite conclusions. 

In the end, let us summarize the main points of the procedure to extract 7 presented 
in this letter. This method brings a direct measurement of the CP violating angles. No 
previous knowledge of BR's or FSI phases is required. No necessity of making complicated 
triangle constructions is needed. The Dalitz analysis deals directly with amplitudes; thus, 
it is linearly sensitive to suppressed decays, as Xc 71 "- Because of all of this, it is natural 
that this method does provide a measurement which demands less statistics than other 
methods. Finally, as one directly gets the angle itself — instead of, e.g. twice its sine or 
its cosine as in other methods — there are no discrete ambiguities. 

The limitations of this method are mainly due to the fact that in order to have a 
complete knowledge of the errors in the extraction of 7 one needs the ratio P/T. One 
hopes that in the near future this ratio will be better known both from the theoretical 
and the experimental sides. As a result, although this method will probably not yield a 
very accurate measurement of 7, it will certainly allow a simple and effective step to have 
a quick constraint in the value of this angle. 

As a by-product of this method to extract 7, we have presented in this letter a general 
procedure to measure CP violating angles in charged three body decays. It is a natural 
and clear method. The procedure is quite general as it applies to any three body decay. 
Moreover, the whole procedure applies for any CP violating phase. For example, using 
this methods with existing data from D meson decays one could obtain upper limits of 
less than 1° for many CP violating angles. This could be used to constrain beyond the 
Standard Model physics. 
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